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A uniform asymptotic theory of the free-space paraxial propagation of coherent 
flattened Gaussian beams is proposed in the limit of nonsmall Fresnel numbers. 
The pivotal role played by the error function in the mathematical description 

"^ ■ of the related wavefield is stressed. 
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1. Introduction 

Several analytical models aimed at describing the propagation of coherent wavefields having 
an "initial" flat-topped profile were proposed in the past. Among them the supergaus- 
sian and the Fermi-Dirac profiles present a very simple mathematical structure but their 
free-space propagation cannot be described without resorting to numerical integration, even 
within paraxial approximation. To overcome such difficulties, in 1994 Gori introduced the 
first model, called flattened Gaussian (FG henceforth) 1], able to describe the paraxial prop- 
agation of coherent flat-topped beams in exact terms [2] . Differently from the supergaussian 
or the Fermi-Dirac, FG profiles are expressed through a finite sum in the following way: 

N 1 

FGU(0 =exp(-aE-^ m ' C 1 ) 

z — ' ml 
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where N denotes the so-called order of the FG profile. Equation (TjQ) was originally derived 
by starting from the identity 

1 = exp(-e 2 ) exp(£ 2 ) , (2) 

and on truncating the Taylor expansion of the second exponential up to N, On further 



rescaling the argument £ by a factor y/N + 1, the FG profile assumes the characteristic flat- 
topped shape which for N = reduces to a Gaussian distribution whereas for N — > oo tends 
to coincide with the characteristic function of the unitary disk. Since their introduction 
FG beams have proved to be a good model for describing the paraxial propagation of flat- 
topped beams, since the initial profile in Eq. (TjQ) can be recast in terms of a finite number 



a. 



of Laguerre Gauss distributions, both of standard and "elegant" kind [3]. Shortly after 
its introduction, it was also proved that the initial distribution in Eq. ([1]) is practically 
undistinguishable from a_ supergaussian one with a suitable and simple adjustement of their 
respective parameters 



[4|. In 2002, Li [5| gave a more rigorous mathematical ground to 
the Gori's trick used in Eq. (J2J), by finding that, for a typical axisymmetric profile, say 
Fn{£)i to display a flat-topped shape it is mandatory to have the first N even ^-derivatives 
to be vanishing at £ = 0. Then, on using such condition, which is certainly satisfied by 
Eq. (JTJ), Li introduced a new class of flat-topped profiles which, differently from the Gori's 
ones, are easily expressable via the superposition of N fundamental Gaussian functions 
having different widths, suitably chosen. Both Gori's and Li's approaches allow the paraxial 
propagation problem of flat-topped beams to be solved without approximations, as the 



resulting fields turns out to be described in terms of finite sums and, at present, Refs. Ill 

n n A 

and 5] have inspired hundreds of papers |6j. At the same time, however, although no 
computational difficulties are encountered in the numerical evaluation of propagated FG 
beams field, the behavior of those sums is hard to grasp, especially from an intuitive point 
of view and this led several authors to think that the FG model has a demanding analytical 
structure. The limiting cases N = and N = oo correspond to initial field distributions 
that are in some respect antithetic as far as their features are concerned: the former is the 
simpler example of tranverse coherent shape- invariant field distribution, whereas in the latter 
the presence of the field discontinuity leads to a complex and mathematically untractable 
(at least from an analytical viewpoint) behavior of the related paraxial free propagation 
features [7| . In between, the FG mathematics, whose exploration is the topic of the present 
work. 

In particular, our aim is to find a uniform [8|, with respect to the transverse radial posi- 
tion, analytical approximation of the free-space propagated field of a iVth-order FG beams. 
The idea for the present work starts just from the limit N = oo which, as said above, 
corresponds to the Fresnel diffraction by a circular hole 3] and from the, apparently not 
so known, asymptotic approximation of the related propagated field found in 1898 by Karl 
Schwarzschild [9J, whose derivation can be found in Sec. 4.3 of [10J], which is a much more 
accessible paper than the Schwarzschild's original one. In particular, within such approxi- 
mation a pivotal role is played, especially for large values of the Fresnel number [7|, by the 
error function which rules the mathematical behaviour of the overall diffracted field. Cer- 
tainly the presence of the error function in the uniform asymptotics of diffracted wavefields 
by hard-edge apertures is not surprising, as pointed out for instance by Stamnes thirty years 
ago [8j, as it makes the field transition across the hole geometric shadow smooth. However, 
although for finite N the FG initial field distribution is continuous, our analysis will show 
that, even for moderately small values of N (of the order of ten), the propagated field is, 
in the limit of nonsmall values of the Fresnel number (roughly greater than two), still well 
represented by an error function. In particular across the initial plane (corresponding to an 
infinte Fresnel number) our analysis provides an interesting and, at least up to our knowl- 
edge, unsuspected direct connection between error function and FG profiles in Eq. ((T]), whose 
mathematical justification will be done a posteriori by using the above quoted Li's flatness 
prescription. 



2. Theoretical Analysis 

We start our analysis from the field distribution of an axisymmetric FG beam of order N 
at the plane z = of a cylindrical reference frame (r, z), given by |l| 

U N {r; 0) = FG^v ( — VnTT) , (3) 

where an unessential amplitude constant has been set to one for simplicity and the symbol 
wo denotes the spot-size of the FG beam. To express the field, say U^{r; z), propagatedjn 
the free space at the transverse plane z > we express the FG profile in Eq. ([I]) as follows 
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where L m (-) denotes the mth-order Laguerre polynomial 11]. The decomposition in Eq. (jlj) 
has a clear physical interpretation once inserted into Eq. ([3]): the field of a iVth-order 
FG beam is thought of as the superposition of elegant Laguerre- Gaussian modes [12J, whose 
functional form is essentially the product of a Laguerre polynomial times a Gaussian function 
with the same complex argument. In this way it_^can be possible to express the propagated 
transverse field distribution f/Ar(r; z) as follows 
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where Np = Wq/\ z denotes the Fresnel number [7(] and where the function Qn\ 
as 

N 



n=0 ^ ' 



is defined 



(6) 



Note, in particular, that in the last equation the variables s and t are non independent, since 



s = (N + l)t. 
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Equations §5§ and ([SD provides the exact expression of the propagated FG beam; moreover, 
the numerical evaluation of the Qn can be efficiently implemented, for large values of N, via 
a simple recursive computational scheme 



To derive a uniform asymptotic approximation of Un{t; z) with respect to r we start from 
the following integral representation of Laguerre polynomials 13, Eq. (5.4.1)]: 

L n {s) = ^[^ |°° d£ exp(-0 r Jo (2 yfit) , (8) 

where J n {-) denotes the nth-order Bessel function of the first kind ll|. On substituting 
from Eq. ([S]) into Eq. §§§ and on taking into account that 
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with Ln 



denoting the generalized Laguerre polynomial of order n and degree a [ll| , the 
following integral representation of the function £?tv(£, s) is obtained: 



exp(-s) G N (t, s) 
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which will be estimated in the limit of non nonsmall values of N. To this end the Laguerre 
polynomial L N is first approximated by using Eq. (8.22.4) of 13j which gives 



cxp 



~o ) 2 Ji 



2y/(N+l)£t 



(11) 



N+1 JV v ^ 2v /(iV + l)^ ' 

that, once substituted into Eq. ( TTUj) . on taking Eq. © into account, and after letting rj 
(N + l)f£, leads to 

exp(-s)Q N (t,s) ~ 
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The integral in Eq. ( 112)) can be evaluated on using formula 2.12.39.1 of [14J 

but straightforward algebra, allows Eq. (J5J) to be recast as follows: 



(13) 
which, after long 
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where the function J{u,p) is defined by 

POO 

J(u, p) = 1 - exp(-wp 2 ) / d£ exp (-£) I (2p y/^u) 

J u 



(15) 



with 7o(") denoting the zeroth-order modified Bessel function of the first kind 11]. 

Differently from the integral representation in Eq. (I23p . the function J~(u,p) is expressed 
in a mathematical form suitable to extract a uniform approximation with respect to p. To 
this end we make use of the, above quoted, Schwarzschild's approach originally employed to 
study the Fresnel diffraction from a circular hole 10J. First of all we replace the modified 
Bessel function by its asymptotic expansion, namely 
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which, once subtituted into Eq. (115j) and after changing the integration variable £ into u £ 2 
gives, for \u\ 3> 1, 
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V£d£exp[- M (e-p) 2 ]. 



(17) 



The last integral can now be estimated by replacing the factor \/£ by (y/£ — y/p) + y/p, so 
that 



J(u,p) ~ 1 - J^ J d£ expi-u (£ - p) 2 ] 
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The first integral can be evaluated exactly by using formula 2.3.15.4 of |14|, which gives 
y^ J™ d£ exp[- M (£ - p) 2 ] = ierfc[v^(l - p)] 

where erfc(-) denotes the complementary error function 11] . As far as the second integral 
in Eq. ( TT8|) is concerned, a first partial integration gives at once 
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while further partial integrations would give rise to an asymptotic expansion in negative 
powers of u whose single terms can easily find by using standard techniques 15] . In Ap- 
pendix |A] it is shown that 

1 r d£ exp[- M (£-p) 2 ] ^ 
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and in the same Appendix it is also given to the interested reader an idea about how to 
derive the complete asymptotic series in such a way that the function J can be written as 
follows: 

J(u,p) = 1 - I- ertc[^(l ~ p)}JT am{p) 
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Equations (122"]) and ( |23l) constitutes the main results of the present Paper and suggest an 
intriguing and, at least up to our knowledge, unespected mathematical connection between 
FG profiles and the error function. 

3. A new analytical form for flat-topped profiles 

First of all we note that in the limit Np — > oo the propagated field Un(v; z) must coincide 
with the initial field Un{t\ 0). From Eq. (fill we have 
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which, once compared with Eq. (J3J), leads to the following asymptotic relation: 
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and where only the leading term in the expansion of Eq. (122]) has been kept. Figure [T] shows 
a visual comparison of FG profiles obtained through the definition in Eq. ([1]) (solid curve) 
and through the asymptotic estimate in Eq. (1251) (dashed curves). 

It is surprising how, except for the case N — 1, the agreement between the two curves 
is remarkably good also for reasonably small values of N (of the order of ten), the two 
curves being practically undistinguishable, at least at a visual level, for N = 30. To give 
an explanation of such behavior it is sufficient to prove that an initial erfc-based profile in 
the r.h.s. of Eq. (|25|) . when is evaluated at the scaled variable £/y/N + 1, satisfies the Li's 
flatness condition which, i.e., 
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In Appendix [B] it is shown that, when the FG profile is approximated via Eq. (|25|) . it turns 
out to be 
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where H n (-) denotes the nth-order Hermite polynomial rjj. Accordingly, we have to prove 
that the r.h.s. of Eq. ( 1271) is negligible for n < N with respect to the values assumed for 
n > N. A possibility is again offered by asymptotics: in Appendix [C] is it shown that, in 
the limit of nonsmall values of n, Eq. (|27|) can be replaced by 
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which, on taking the absolute value of both members and on replacing the modulus of the 
sinusoidal function by its unitary upper bound, leads to the inequality 
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Figure [2]^a) shows a visual comparison between the behaviors of the modulus of the 2nth- 
order ^-derivative in Eq. ( 127)) (open circles), together with the asymptotic upper bound 
given in Eq. (12"9~)) (solid curve) for a typical value of the FG order iV (in this case iV = 30). 
In order to better appreciate, from a quantitative point of view, the agreement between 
the two behaviors, the same plot is drawn in Fig. |2J(b) but on a vertical logarithmic scale. 
Inequality (129)) is the key to grasp why the erfc-based approximation in Eq. (125)) satisfies 
the Li's flatten condition: it predicts a "power times factorial" law which is known to 
displays the characteristic treshold-like behavior of Fig. [2j(a) [l6(. In particular, the terms of 
the sequence corresponding to the r.h.s. of Eq. (12"9")) turns out to be exponentially small for 
n < N, whereas they grow according to a factorial divergence law for n > N. In Appendix ID1 
a nonrigorous proof that the treshold is around n ~ N is given to the interested reader. 

4. Some numerical results on free-space propagation 

The present section is aimed at exploring, from a numerical point of view, the free-space 
propagation of FG beams by presenting a comparison between the exact values, obtained 
through Eqs. (jSJ) and (Q, and those provided by the uniform approximation in Eqs. (fl4j) . (1221) . 
and (123)) for a FG profile of order N = 30. In Fig. [3] the amplitude (a) and the phase (b) 
of the FG beam propagated at the Fresnel number Np = 10 are shown as functions of the 
normalized transverse radial distance r/wo (the dotted curve is the initial FG profile). Open 
circles are the exact values provided by Eq. §J5§, whereas the solid curve is the asymptotic 
estimate obtained by keeping only the leading term in the asymptotics expansion of Eq. ( 1221 ). 
The dotted curve represents the (normalized) initial FG profile. As we can see, the agree- 
ment is really good. On reducing the Fresnel number to Np = 5 the amplitude and phase 
distributions are those reported in Fig. H^a) and Fig. H](b), respectively. Also in this case it 
is possible to appreciate a reasonably acceptable agreement, especially far from the beam 
axis. On further reducing the Np the agreement continues to dedragates: the situation at 
Np = 2 is depicted in Fig. [5] where it is clear that the sole leading erfc-based asymptotics is 
no longer able to guarantee an acceptable, even at a visual level, agreement for the trans- 
verse beam amplitude. To try to improve such an agreement we add to the leading term of 
the asymptotics expansions in Eq. ( 122)) that corresponding to bo, thus obtaining the plots of 
Fig. El which now displays a much better agreement, especially close to the beam axis. 
To conclude the present section we want to present a simulation concerning with the 



case Np = 1, which is customarily assumed as the boundary between the near-zone and 
the far-zone [2J. Figure [7] shows the amplitude (a) and phase (b) distributions obtained by 
keeping only the leading erfc-based term (dotted curve), by including the term b (dashed 
curve), and by adding also a\ and a<i (solid curve). It must be appreciated that, although 
our analysis has been formally developed in the limit Np ^> 1, the asymptotic series found 
above are able, with a modest increase of the computational effort, to provide a reasonably 
good description of the propagated field also beyond the analysis limits. 

5. Conclusions 

Still nowadays the supergaussian seems to be, especially among experimentalists, the most 



known and used model for dealing with flat-top beams [17H21|. although it is well known 
that their propagation features are untractable in analytical terms. On the other hand, the 
"exact" models introduced by Gori and later by Li suffered by a (only apparent) formal 
complication of the propagated field expression which led several authors to still prefer 
the older and obsolete models. In the present paper we have shown that some interesting 
general features of the free-space propagation of flattened Gaussian beams can be grasped 
via uniform asymptotics. On using the FG beam representation in terms of elegant Laguerre- 
Gauss modes together with the Schwarzschild's approach to derive uniform approximations 
of integrals we have found that the propagated field is well approximated, for large Fresnel 
numbers, by a simple error function. In particular, starting from the initial plane, the 
functional form of the propagated field remains the same and the effect of the propagation is 
taken into account by a simple (complex) scaling of the error function argument. On further 
reducing the Fresnel number this sort of "functional invariance" is progressively lost and 
it is necessary to include higher-order terms of the asympotic series to achieve acceptable 
agreements between the exact and the approximated field. Since it has been proved in the 
past the practical indistinguishibility of FG and supergaussian profiles, the results obtained 
in the present paper apply, upon evaluation of the related parameters, also to the latter. 
Moreover, a byproduct of our analysis is a new, at least up to our knowledge, analytical 
expression for flat-topped profiles, whose related flatness conditions have been quantitatively 
verified according to the Li's prescriptions. 

Finally, while the role of the error function in the mathematical treatment of diffraction 
by hard-edge apertures is well known, our results seem to suggest its involvement also 
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when dealing with smoother initial field distributions; moreover, since the representation 
of FG beams via elegant Laguerre- Gauss modes is not limited to the case of the free-space 
propagation it should be possible, in principle, to extend the analysis presented here to the 
propagation problem through a typical paraxial ABCD optical system. 

Appendix A: Proof of Eq. ( 12T1) 

The Schwarzschild's trick can be thought of as the first step of a systematic procedure to 



extract uniform approximations of integrals, as pointed out for instance by Temme 15] . To 
show this consider first the evaluation of integrals of the type 

d*/(0 exp[-w(£ - pf] , (Al) 
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where /(•) denotes a function sufficiently regular in the interval [1, oo]. To obtain a uniform 

approximation with respect to p we have to replace the factor /(£) with f(p) + [/(£) — /(/>)] 

in such a way that Eq. (IA10 becomes 
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The first integral can be evaluated again via Eq. (Tl9|) . while the second integral can be 
rearranged for a further partial integration, so that 
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where the last integral gives a contribution of order 1/u 2 . To prove Eq. (12TJ) it is sufficient 
to substitute into Eq. flA3j) the function /(£) given by 

m = vmrw- (A4) 



thus obtaining 
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where the last term can be further simplified on taking into account that 
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so that, after rearranging, Eq. (1A5J) becomes 
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which, on neglecting the last integral, coincides with Eq. (J2T 
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Appendix B: Proof of Eq. ([2] 

On substituting from Eq. (125]) into Eq. (126]) and on letting ( = £,/\/N + 1, we have 
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and, on using formula 1.5.1.1 of 22], Eq. ( 127|) follows. 
Appendix C: Proof of Eq. (|28|) 

We start from the following asymptotics of Hermite polynomials [see formula 8.22.7 of [131]]: 
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valid for nonsmall values of n. On using Stirling's formula to estimate the factorials, Eq. ljClj) 
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which, once substituted into Eq. ( 12"T|) . after some algebra leads to Eq. (128]) . 
Appendix D: Proving that the threshold is at n ~ iV 

We could conventionally set the treshold as the value of n at which the modulus of the r.h.s. 
of Eq. ( 129]) is of the order of the unity. Taking the logarithm this implies that, neglecting 
the factor l/(2ir), 



logT(n) — nlog 



JV + 1 



(JV + 1) ~ 0, n>l 
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where T(-) denotes the gamma function Uj. For large values of n the following asymptotics 
holds: 

logr(ri) ~ nlogn — n , n^$> 1 , (D2) 
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which, after substituted into Eq. ( IDlj) and after taking into account that log 8 ~ 2, gives 

nlogn +n - n\og(N+ 1) - (N + 1) ~ 0, n > 1 , (D3) 

which, in the limit of nonsmall iV's, gives n ~ TV. 
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List of figure captions 

Fig. 1. Flattened Gaussian profiles FGjv(£) evaluated, for different values of N, through 
Eq. ([I]) (solid curves) and through the asymptotic estimate in Eq. (^5§ (dashed curves). 

Fig. 2. (a): behavior of the modulus of the 2nth-order ^-derivative in Eq. (|27[) (open circles), 
together with the asymptotic estimate in Eq. (I2"9j) (solid curve). N = 30. (b): the same plot 
as in figure (a), but on a vertical logarithmic scale. 

Fig. 3. Behavior of the amplitude of a FG beam of order N = 30 propagated at the Fresnel 
number Np = 10. Open circles: exact values provided by Eq. (jol) . Solid curve: erfc- 
based asymptotic estimate by keeping only the leading term in the asymptotics expansion 
in Eq. (I22p . The dotted curve represents the (normalized) initial FG profile. The phase 
distribution is wrapped. 

Fig. 4. The same as in Fig. [3j but for Np = 5. 

Fig. 5. The same as in Fig. [3j but for Np = 2. 

Fig. 6. The same as in Fig. but now including the term 6 i n the asympotic expansion of 
Eq. ([22]). 



Fig. 7. Amplitude (a) and phase (b) distributions obtained by keeping only the leading 
erfc-based term (dotted curve), by including the term bo (dashed curve), and by adding also 
ai and ci2 (solid curve). Np = 1. 
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